Penetration depth, symmetry breaking, and gap nodes 
in superconducting PrOs 4 Sb 12 
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Transverse-field muon spin relaxation rates in single crystals of the heavy-fermion superconductor 
PrOs4Sbi2 (T c — 1.85 K) are nearly constant in the vortex state for temperatures below ~0.5T C . 
This suggests that the superconducting penetration depth A(T) is temperature-independent at low 
temperatures, consistent with a gapped quasiparticle excitation spectrum. In contrast, radiofre- 
quency measurements yield a stronger temperature dependence of A(T), indicative of point nodes 
in the gap. A similar discrepancy exists in superconducting Sr2Ru04 which, like PrOs4Sbi2, breaks 
time-reversal symmetry (TRS) below T c , but not in a number of non-TRS-breaking superconductors. 

PACS numbers: 71.27. +a, 74.70.Tx, 74.25.Nf, 75.30.Mb, 76.75.+i 



Phase transitions are always symmetry breaking; the 
symmetry of the ordered state is less than the full sym- 
metry of the Bamiltonian. In all superconductors gauge 
symmetry is broken and the ground-state wave function 
acquires a definite phase. Additional symmetries are si- 
multaneously broken in so-called unconventional super- 
conductors. Examples include the point-group symmetry 
of the lattice, spin rotation symmetry, and time-reversal 
symmetry (TRS) (HQ- TRS is always broken in mag- 
netic transitions but rarely in superconductors. Prob- 
ably the most direct evidence for TRS breaking is the 
onset of a spontaneous local magnetic field in zero ap- 
plied field below the superconducting transition temper- 
ature T c ; the local field reverses its direction under time 
reversal and thus signals broken TRS. 

Muon spin rotation/relaxation (/^SR) experiments 
have proved invaluable in characterizing the supercon- 
ducting state Q . In the /iSR technique spin-polarized 
positive muons (/i + ) are implanted in the sample and 
precess in the local magnetic field. This precession is 
detected using the asymmetry of the beta decay (the 
decay positron is emitted preferentially in the direction of 
the /x + spin) . The distribution of fi + precession frequen- 
cies directly reflects the distribution of magnetic fields 
in the sample. Reproducible ^iSR evidence for a TRS- 
breaking local field has been found in only three super- 
conductors: the heavy-fermion system (U,Th)Bei3 H, 
the transition-metal oxide S^RuO^ 01, and the Pr-based 
heavy-fermion compound PrOs4Sbi2 |2j- 

TRS breaking in a superconductor may (but need not) 
be accompanied by breaking of other symmetries 0. 



The most accessible consequence of additional symmetry 
breaking is the presence of point or line nodes A(k) = 
in the superconducting energy gap function A(k) on the 
Fermi surface. Thermal excitations of low-lying states 
are strongly modified by gap nodes, and these excita- 
tions govern superconducting properties at low tempera- 
tures. The temperature dependence of the superconduct- 
ing penetration depth A, in particular, depends on details 
of the gap structure 0, Q : for a nodeless superconductor 
AA(T) = A(T) - A(0) oc e" A / T , which is small at low 
temperatures, whereas with gap nodes AA oc T n , n < 2, 
for at least some field orientations. 

PrOs4Sb 12 is a superconductor (T c = 1.85 K Q) with 
a number of extraordinary properties [jj . It is the only 
known Pr-based heavy-fermion superconductor, the Pr 3+ 
ground state is nonmagnetic, a novel ordered phase ap- 
pears at high fields and low temperatures, there are mul- 
tiple superconducting phases, and TRS is broken in the 
superconducting state Q. Previous /iSR measurements 
of A in the vortex state of a powdered sample [lfj found 
evidence for a BCS-like activated dependence at low tem- 
peratures, suggesting the absence of gap nodes. But ra- 
diofrequency (rf) inductive measurements of the surface 
penetration depth indicate point nodes in the gap, 
in disagreement with the fj.SK results. 

This Letter reports new /iSR experiments on ori- 
ented PrOs4Sbi2 crystals, and compares /iSR and sur- 
face penetration-depth measurements in PrOs4Sbi2 and 
a number of other superconductors. The discrepancy 
between these measurements noted above is found in 
TRS-breaking superconductors but not otherwise, and 
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is therefore correlated with TRS breaking. 

In the vortex lattice of a type-II superconductor each 
vortex possesses a normal-state-like core surrounded by a 
shielding supercurrent. These supercurrents give rise to 
an inhomogcneous magnetic field that is periodic in the 
vortex lattice. The vortex-lattice field distribution P V (B) 
is related to A. For a fixed density of vortices the rms 
width SB of P V (B) decreases with increasing A; in a Lon- 
don superconductor (A much longer than the coherence 
length £) SB oc 1/A 2 when the applied field H satisfies 
H cl < H < H c2 [13 Transverse-field /zSR (TF-^SR) 
experiments (field applied transverse to the spin) are 
sensitive probes of this field distribution Q . 

Time-differential TF-/iSR experiments were carried 
out at the M15 channel at TRIUMF, Vancouver, Canada, 
on a mosaic of oriented PrOs4Sbi2 crystals. The crystals 
were mounted on a thin GaAs backing, which rapidly 
depolarizes muons in transverse field and minimizes any 
spurious signal from muons that do not stop in the 
sample. /iSR asymmetry data were taken for tem- 
peratures in the range 0.02-2.5 K and fi^H between 
10 mT and 100 mT applied parallel to the (100) axes 
of the crystals. The data were fitted with the functional 
form G(t) cos(ujt + (j)) 1 where the frequency lo and phase <j) 
describe the average fj, + precession and the relaxation 
function G(t) describes the loss of phase coherence due 
to the distribution of precession frequencies 0. The 
relaxation rate associated with G(t) is a measure of the 
width of this distribution and hence of SB [j| . 

Neither of the commonly-used exponential or Gaus- 
sian functional forms accurately fit the asymmetry data 
in the normal state, due to // + coupling to both nuclear 
and Pr 3+ spins as discussed below. Data from both the 
normal and the superconducting states are well fit, how- 
ever, by either of two slightly more complex functional 
forms: the "power exponential" 
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G(<) = cxp[-(At)*], 



and the damped Gaussian 
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These functions are both phenomenological and have no 
theoretical motivation. We shall see that the supercon- 
ducting-state properties obtained from these fits are sim- 
ilar for both functions, which indicates insensitivity to 
details of the fitting function and justifies a posteriori 
these otherwise arbitrary choices. 

Figure H^a) shows the temperature dependence of the 
relaxation rate A and exponent K from power-exponent- 
ial fits for ii H = 10 mT and 100 mT. The relax- 
ation is intermediate between exponential and Gaussian 
(1 < K < 2). The temperature dependence of K below 
T c indicates that the shape of G(t) is changing with tem- 
perature, but for temperatures lower than about 1 K A 
and K are nearly constant for (IqH = 10 mT. 
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FIG. 1: Relaxation rates from TF-/iSR asymmetry data in 
PrOs4Sbi2 for applied fields fioH = 10 mT (circles) and 
100 mT (triangles), (a) Relaxation rates A (open symbols) 
and exponents K (filled symbols) from power-exponential fits 
[Eq. JJJ], Note that A and K have different dimensions; they 
are plotted on the same graph only for convenience, (b) Gaus- 
sian relaxation rates a from damped Gaussian fits [Eq. JSJ]. 
The exponential rates have been fixed to the normal-state 
values [W(10 mT) = 0.095 /as -1 , W(100 mT) = 0.428 ^s" 1 ]. 



Figure ^b) gives the temperature dependence of the 
Gaussian rate a from damped-Gaussian fits. Below T c 
the exponential rate W [Eq. was held fixed at the 
normal-state value for each field, so that the tempera- 
ture dependence of the Gaussian rate a reflects the effect 
of the superconducting state. Some such procedure is ne- 
cessitated by the strong statistical correlation between W 
and a in Eq. (J2J; the time constant and the shape of the 
function are influenced by both of these parameters, so 
that correlations between them can result from small sys- 
tematic errors. The principal justification for this ad hoc 
fixing of W is the insensitivity of the superconducting- 
state results to details of the fitting function noted above. 
Statistical correlation is less of a problem in power expo- 
nential fits, where A and K control the rate and shape, 
respectively, and are not as strongly correlated. 

For damped-Gaussian fits the normal-state Gaussian 
rate is independent of field [cf. Fig. IHb)] and of mag- 
nitude consistent with nuclear dipolar broadening. The 
exponential rate increases with field and is due to in- 
homogeneity in the Pr 3+ magnetization. Determination 
of the vortex-state field distribution width requires cor- 
rection for the normal-state relaxation. For the damped 
Gaussian fits we take the superconducting-state Gaus- 
sian rate a s to be given by c 2 = a 2 — cr 2 , where a n is the 



normal-state rate 

3]- 

For the power exponential fits the 
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proper correction procedure is somewhat less clear. We 
have chosen to use the relation A^ = A K — A„ , which 
interpolates smoothly between the exponential (K = 1) 
and Gaussian (K = 2) limits. The results are sensitive 
to this choice only for /irj-ff = 100 mT, where A n is sig- 
nificant [cf. Fig. I^a)]. 

Figure [21 gives the temperature dependence of the 
corrected superconducting-state fi + relaxation rates for 
HoH = 10, 20, and 100 mT. It can be seen that the 
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FIG. 2: Temperature dependence of superconducting-state re- 
laxation rates in PrOs4Sbi2, corrected for normal-state relax- 
ation (see text), (a) power exponential rates A s . (b) damped 
Gaussian rates a 3 . Circles: applied field noH = 10 mT. 
Squares: fi H = 20 mT. Triangles: ^jl H = 100 mT. The 
curves are guides to the eye. 

qualitative behavior of the rates is remarkably indepen- 
dent of the fit function. At 10 and 20 mT both A s 
and a s are nearly temperature- independent below ~1 K. 
At the lowest temperatures the rate is field-independent 
to within a few percent between 10 mT and 100 mT. 
In an isotropic superconductor such as cubic PrOs4Sbi2 
vortex-lattice disorder is expected to increase the low- 
field rate; increasing field (increasing vortex density) then 
decreases the rate as intervortex interactions stabilize the 
lattice 0, 0] . Thus the field independence of the rate 
indicates a substantially ordered vortex lattice, in which 
case the temperature dependence of the rate is controlled 
solely by the temperature dependence of the field distri- 
bution P V (B). We can also conclude that the field de- 
pendence expected as H — > H cl 15j plays no role. 
The expression 



5B 2 (T) = 0.00371 $o A " 4 ( T ) ■ 



(3) 



where $o is the flux quantum, relates the second mo- 
ment SB 2 of P V (B) to A for a triangular vortex lattice 



in the London limit. The second moment of the cor- 
responding /i + frequency distribution is 5lo 2 = j 2 SB 2 , 
where 7^ is the /i + gyromagnetic ratio. Then the /xSR 
estimate A m sr of the penetration depth from Eq. J5J is 



Ksr (Mm) = 0.328/vAW (^s- 1 ) . 



(4) 



Now the rms width a s of the best-fit Gaussian is not 
necessarily 8uj, and replacement of 5tu in Eq. Q by a s 
is not completely justified. Nevertheless a s should scale 
with 5u>, and within its range of validity Eq. Q should 
give the correct temperature dependence of A^sr- This 
is because under these circumstances effects of nonzero 
£ are restricted to the high- field tail of P V (B), which 
is not heavily weighted in a Gaussian fit (cf. Fig. 1 of 
Ref. 0). P rOs4Sbi2 is a strongly type-II superconduc- 
tor (Ginzburg-Landau k = A/£ « 30 H0), and this 
picture should be applicable. 

Figure Ofa) compares AA surf (r) = A surf (T) - A surf (0), 
obtained from rf inductance measurements in the Meiss- 
ner state 0,0], with AA M sR(r) obtained using /i+ re- 
laxation rates for /i 7f = 10 mT (Fig. EJ) in Eq. J3J. 
At low temperatures 0] the difference AA sur f(T) — 



0.02 
0.01 
0.00 

0.04 
0.02 
0.00 
0.1 
0.0 



E 



< 



< 




(d)YNLBC 



♦ T 



(b) Sr 2 Ru0 4 



(e)YBa 2 Cu 3 O eM 



-0.1 



(c) CeColn 6 



0.00 • 



-0.05 



WLa^Sr^CuO, 



0.0 0.2 0.4 0.6 
TIT 



0.0 0.2 0.4 0.6 
777 



FIG. 3: Dependence of difference AA sur f — AA m sr (see text) 
on reduced temperature T/T c in six superconductors, (a) 
PrOs4Sbi2. Triangles: A^sr from A s [Fig. Ufa)]. Circles: 
A mSR from a 3 [Fig. Hb)]. (b) Sr 2 Ru0 4 . (c) CeCoIn 5 . (d) 
YNi 2 B 2 C. (e) YBa 2 Cu 3 6 .95. (f) Lai. 85 Sro.i 5 Cu04. See text 
for references. 

AA m sr(7^) increases markedly with increasing tempera- 
ture; this is the discrepancy between the measurements 
noted above. It is the same whether A S (T) [Fig. [^a)] 
or ct s (T) [Fig. |2Jb)] is used in Eq. |0J; the compari- 
son does not depend on the choice of fitting function. 
Figure Hb) gives AA surf (T) - AA AlSR (T) for the TRS- 
breaking superconductor Sr2Ru04, using data from the 
literature 0, 0|. Again there is a discrepancy (as 
noted previously |2l|). which is very similar to that 
in PrOs^Sbi2. Small-angle neutron diffraction experi- 
ments 22] also found a temperature-independent vortex- 
lattice field distribution in Sr2Ru04 at low temperatures. 
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Figures Etc) 
from literature 



-(f) give 
data for 



AA surf (T) - AA mSR (T) 
the heavy-fermion com- 




pound CeCoIn 5 [23, [2J|, the borocarbide YNi 2 B 2 C 
l2q ]. and the high-T c cuprates YBa 2 Cu 3 6 .g5 
and Lai.85Sro.i5Cu04 [23, H(J. None of these super- 
conductors exhibit TRS breaking, and none exhibit 
the temperature dependence of AA sur f — AA^sr seen 
in Figs. Eta) and (b). There is a significant difference 
between AA sur f(T) and AA m sr(T) on ^V f or ^ e TRS- 
breaking superconductors PrOs^Sb^ and Sr2RuO^. This 
is in spite of a wide range of experimental methods 
and analysis techniques: in Figs. 3(a), (b), (c), and (f) 
AA sur f was obtained from rf measurements, in (d) from 
magnetization, and in (e) from microwave impedance; in 
Figs. 3(a) and (c) A m sr was obtained from Eq. Q), in 
(b), (e), and (f) from fits to the expected vortex-state 
field distribution, and in (d) from a generalization of 
Eq. <gj) to low and high fields |!|. 

The origin of this discrepancy and its relation to TRS 
breaking are not clear. Low-field low-temperature phase 
transitions between superconducting states have been re- 
ported in both PrOs4Sbi2 and Sr 2 Ru04 and 
may be involved in the discrepancy. The TRS-breaking 
state may couple to rf or microwave fields, necessitating 
a revised interpretation of the surface measurements. A 
mechanism of this sort, in which subgap chiral surface 
states affect the surface penetration depth, has been pro- 
posed for Sr 2 Ru04 33]. A T 2 power law is found for 
surface measurements even though the bulk energy spec- 
trum is gapped. The theory requires A » £, however, 
and thus seems inapplicable to PrOs4Sb i2 . It has also 
been noted [34| that surface scattering breaks pairs in an 
odd-parity superconductor. To our knowledge the surface 
penetration depth has not been calculated taking this ef- 
fect into account, but pair breaking would decrease the 
gap and therefore increase the temperature dependence 
of A. The discrepancy might be related to a breakdown 
of the relation 5B(T) oc 1/A 2 due to nonlinear/nonlocal 
effects UE3j_P r to tne spontaneous magnetic field in the 
vortex state pf (although the measured field is not large 
enough to have a significant direct effect on SB). 

Sr 2 RuC>4 and PrOs4Sb i2 are both TRS-breaking su- 
perconductors but are otherwise very different. Tetrago- 
nal Sr 2 RuC>4 is an anisotropic transition-metal-oxide su- 
perconductor that is weakly type-II for H || c {n a b = 
2.3 whereas cubic PrOs4Sb i2 is an isotropic 

strongly type-II heavy-fermion superconductor. The fact 
that a similar discrepancy between ^tSR and surface mea- 
surements of A(T) is found in such different materials, 
but not in a variety of non-TRS-breaking superconduc- 
tors, strongly suggests that TRS breaking is involved. 
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